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Abstract 

We study the non-equilibrium relaxation of the spherical spin-glass 
model with p-spin interactions in the N — > oo limit. We analytically 
solve the asymptotics of the magnetization and the correlation and re- 
sponse functions for long but finite times. Even in the thermodynamic 
limit the system exhibits 'weak' (as well as 'true') ergodicity breaking 
and aging effects. We determine a functional Parisi-like order param- 
eter Pd(q) which plays a similar role for the dynamics to that played 
by the usual function for the statics. 
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Spin-glass dynamics has been a subject of continuous interest in the past 
years. Experimentally, spin-glass dynamics below the critical temperature is 
characterized by aging effects and very slow relaxations []T|. In long range 
mean-field models one knows that the phase space is broken into ergodic com- 
ponents ||. Sompolinsky || described a dynamics for these models allowing 
for barrier penetration in very long times (diverging as N — > oo). 

In realistic systems, on the one hand mean-field is not exact and on the 
other hand one cannot perform an experiment in infinite times, and one 
actually sees at most 'weak' ergodicity breaking. 

Bouchaud has proposed a phenomenological scenario with both 'true' 
and 'weak' ergodicity breaking [|J. The question then arises as to if and 
how simple long-range microscopic systems (for which mean- field is exact) 
can model these phenomena. To the best of our knowledge, an analytic 
description is lacking. 

The main purpose of this paper is to show, in a very simple mean-field 
model, the asymptotics of which we solve analytically, that this is indeed 
so; in the thermodynamic limit 'true' and 'weak' ergodicity breaking coexist, 
and in a sense complement. To this end we solve the dynamics of the p-spin 
spherical model (p > 2) first introduced in ref. I|, setting N — > oo from the 
outset, starting from a given configuration, for long (but not diverging with 
N) times. 

It should be stressed that this is a different physical situation from the 
Sompolinsky dynamics, which was analysed in ref.||. We do not have here 
any time-scale dependent on iV (or any other 'regularization' parameter): 
for the two-time (t, t') functions the scale that naturally arises is t'/t. Sur- 
prisingly, one can establish formal contact with Sompolinsky's equations by 
defining a variable r = log(t'/t) which plays the role of the 'time' there (this 
will be further explained in a separate work in the context of the Sherrington- 
Kirkpatrick (SK) model [§). 

Denoting the two-spin correlation and the linear response to a magnetic 
field 



N i N 



d(s t (t)) 



our main results £1X6 cLS follows: 

i. For any waiting time t w there exists a sufficiently large t such that 
C(t + t w , t w ) tends to zero. 
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ii. After a (large) t w , the decay rate of the correlation function has a 
short transient after which it is inversely proportional to t w (an aging-like 
effect). 

in. For t large, the magnetization falls to zero as t~ u . 

iv. In addition to a strong short-term memory, the system possesses a 
weak, long-term memory. 

We expect other models, such as the Potts glass (for more than three 
components) and the p-spin Ising model (for not too low temperatures) to 
have a similar dynamics to the one presented here. The SK model instead 
has a rather different behaviour M . 

The spherical p-spin glass model is defined by the Hamiltonian 

N I N 

H = — E Ji 1 ...i p Si 1 ■ ■ ■ Si p + _ x E ^ii...ir S n • • • S ir • (1) 
ii<...<i p ii<...<i r 

The spin variables verify the spherical constraint X^i 5 ?^) — N. The in- 
teraction strenghts are independent random variables with a Gaussian dis- 
tribution with zero mean and variance (J^...^) 2 = p\/(2N p ~ 1 ). The overline 
stands for the average over the couplings. Additional source terms (/i^...^ 
time-independent) have been included; if r = 1 the usual coupling to a mag- 
netic field hi is recovered. 

The relaxational dynamics is given by the Langevin equation 

r 1 d tSt (t) = -Pj^z - + • ( 2 ) 

T determines the time scale and will be henceforth set to one. The second 
term in the rhs enforces the spherical constraint while is a Gaussian 
white noise with zero mean and variance 2. The mean over the thermal noise 
is hereafter represented by ( ■ ). As will be shown below, the dynamical 
equations plus the spherical constraint impose z{t) = (1 —p(3£(t)) with £{t) 
the energy per spin. We choose as initial configuration Sj(0) = 1 Vz, though 
any other choice is equivalent. 

The mean-field sample-averaged dynamics for iV — > oo is entirely de- 
scribed by the evolution of the two-time correlation and the linear response 
functions. The dynamical equations for them can be obtained from eq.(0) 
through standard functional methods (see eg. ref.||) 

^Ml = -(l- pi 3S(t))C(t,t') + 2G(t',t) 
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+ nf t dt" C p -\t,t')G{t\t") 
Jo 

+ im {p - 1) f dt" G{t, t") C p ~ 2 {t, t") C{t", t') , (3) 
Jo 

-(l-pf3£(t))G(t,t') + 5(t-t') 

+ fi(p-l) f dt" G(t, t") C p ~ 2 (t, t") G(t", , (4) 



dG(t,t') 
dt~ 



with \i = p/3 2 /2. These equations hold for all times t and t'. At equal times 
C(t,t) = 1, linv^ t - G(t, t') = 1 and lim t /^ t ± d t C(t, t' ) = ±1. £{t) can be 
identified as the energy per spin multiplying eq.(^) by Sj(t'), averaging over 
the noise and the couplings and taking the limit t' — > t. Furthermore, with 
the definition 



E d< Sl ^ h -- Sdt) > \ h ^=rfdfC-\t,f)G{t,f) 

N hZ^<i r d K-i r JO 

(5) 

eq.flSfl implies £{t) = {(3/2)I p (t). 

In order to make the solution to these equations intelligible, we first briefly 
describe the structure of the TAP free-energy landscape 0, though we shall 
never use the TAP results in the dynamic treatment. The TAP free-energy 
can be written in terms of §i = rrii/ ^fq and q = (1/-/V) m i where rrii are 
the magnetizations 

/tap = <fi Eo{h, ...J N )-^s ln(l - q) - | [(p - l)<f - W^ 1 + 1 

£o(£i, • • • , Sat) denotes the zero temperature energy of a configuration {.§;}. 
The free-energy landscape in the 'angular' variables s, is unaltered by tem- 
perature apart from a stretching proportional to q p ^ 2 . 

The 'angular variable' saddle point equations are supplemented by the 
condition of minimization with respect to q 

q ^- 1 (i-q) = ^- 1 [-£o + (£ 2 -£l) 1/2 } (6) 

where £$ c = — {2{p — l)/p) l ^ 2 \ the largest root for q corresponds to the 
minimum. 
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The paramagnetic solution q = exists for all temperatures. Each zero 
temperature saddle point {§i} with energy £ determines the temperature T 
saddle point (m 8 = Sj^/g, q) where q is obtained from eq.(^). Thus, all TAP 
saddle points are labelled by their associated zero temperature energy £ . 
Moreover, it is easy to see that their ordering in free-energy does not change 
with temperature (i.e. there is no 'chaoticity' with respect to temperature 
in this model). 

Above the threshold value £q c for £ Q (corresponding to a threshold for 
/(/?)) eq.(||) has no solutions. Since it will turn out that the dynamics is 
dominated by this threshold level, it is useful to describe it in more detail. 

Using standard methods ]l(J one finds that the typical spectrum of the 
free-energy Hessian in a local minimum corresponds to a 'shifted' semicircle 
law, with the lowest eigenvalue \ m i n given (in terms of the parameters of the 
minimum) by 

Amin = pq 2 (£oc ~ £q) ■ 

Hence, for sub-threshold free-energies we have well-defined minima with no 
'zero-modes' separated by O(N) barriers. In particular, this was shown 
within the replica approach for the lowest minima that dominate the Gibbs- 
measure ||. Exponential decays would be expected within them; however, 
those low-lying states are quite irrelevant for the non-equilibrium dynamics 
of this model. The gap X min drops to zero at the threshold, and around that 
value the barriers drop from O(N) below to zero above. 

The parameter q and the TAP energy at the threshold are given by 
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/^tk 2 (1 - to) 2 , (7) 



p — 1 

£th 
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i - id - (p ' 2)(1 ' te) ) 

Qth 



(8) 



Let us now turn to the solution of the dynamical equations. Since we 
are interested in the non-equilibrium dynamics we solve them with the only 
assumption of causality. We take t > t' for definiteness and we focus on the 
low temperature phase. The system can be solved numerically step by 



step in a manner reminiscent of ref. |TT[. The numerical solution suggests the 



following scenario for the asymptotic regime t » 1 which we later confirm 
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analytically. The time axis if is divided in three distinct zones with different 
behaviours. 

i. If t' is close to t but (t — t')/t — > asymptotically, time homogeneity 
and the fluctuation-dissipation theorem (FDT) hold; i.e. Gfdt(£ — t') = 
—Q(t — t') dtC-pmit — t'). For large values of (t — f) (but still small compared 
to t) Cfdt(£ — f) tends to a value q and Gfdt(^ — f) tends to zero. 

it. If if is such that (t — if) ft ~ 0(1), the relevant (adimensional) inde- 
pendent variable turns out to be A = t'/t (0 < A < 1). In this sector the 
correlation and rescaled response functions depend on A as C(t,t') = qC{\) 
and tG(t,t') = G(\). Since q is the limiting value of C(t,t') in the previous 
regime, C(l) = 1. 

in. Finite times t' correspond to A = in rescaled variables. In particular, 
for t' = we have the magnetization m(t) = C(t,0). 

We now proceed to solve the resulting equations within this asymptotic 
scenario. If t' is such that the system is in the FDT regime eq.@ yields 

(^- + l)C Fm (t)+(ii+pf3£ 00 ) (l-CWW) = A* Jl dt" C^Ut-t") 

with the asymptotic energy (nb. 'oo' is understood as a limit taken after 
N — > oo) given by 
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'1 - q p ) + pq?- 1 I' d\" g(X")C p - l (X'' 
Jo 

The correlation decays to a value q determined by 



(9) 



1 - pfiSn + y. (1 - g^ 1 ) = . (10) 

This equation appears in the dynamics a la Sompolinsky of this model |J. 
The solution for q as well as the decay law requires solving the coupled system 



(|9T)-(|T0"1) which involves the previous history through the A-integration. 

We now consider the regime < t'/t < 1. The dynamic equations for this 
range of times reduce to two coupled equations for C(A) and £7(A) in which, 
consistently, all times enter only through A: 

= ^(A)[-(l-g)- 1 + / i(l-g)(p-l)g^ 2 ^- 2 (A) 
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C(A) [-(l-g)- 1 + / ,(l-g)g ? '- 2 C p - 2 (A); 

+H (p - 1) g p ~ 2 jT 1 dA" 0(A") C p - 2 (A") C f( A)^(a"-a) 



(11) 



(12) 



Eq.([TT|) in A = 1 admits the solution G(l) = which implies Q(X) = and 
this is the high temperature asymptotics. In the low temperature phase a 
non-trivial G(X) is possible provided the first square bracket in (|TT| ) evaluated 
in A = 1 is zero; this fixes the value q. From eq.([TT|) it also follows GiX) = 
xqC'(l) (prime denotes derivative with respect to A) x = (p — 2)(1 —q)/q. It 
is now easy to see that the system (p|) - (p^ ) with G(X) = xqC'(X) simplifies 
to a single equation. With this ansatz the system of equations has the unique 
family of (exact) solutions 



C(A) 



A 1 



C(M') 



In order to determine z/ (0 < v < 1) a careful matching between this solution 
and the ones associated with other sectors has to be made. 
All the integrals fl5|) in the large t limit become 



1 -q r {l -x) 



(13) 



In particular, E^ = {(5/2)1^. We have now everything that is required to 
solve the FDT relaxation, eqs.(3)-([l0|), which for this value of E^ imply a 
power law decay for this regime ||. Interestingly enough, the expressions 
just derived for the energy E^ and q coincide with eqs. ([?])- (||). Therefore, we 
have learned that the long time dynamics takes place in the threshold of the 
TAP free-energy. 

Finally, we consider the finite t' regime. We already know that for large 
t correlations relax to zero; we now study the asymptotics. Inserting the 
behaviour C(t,t') ~ t~ a c(t') in eqs. (§)-(§) and using the previous results, we 
find a = is, i.e. the exponent for t is the same as in the previous regime. 

The numerical solutions show that the asymptotic regime is well estab- 
lished already for (adimensional) times t ~ 100. 
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We have hence the following picture. For low temperature and increasing 
times the system explores deeper and deeper traps, the permanence time in 
a trap being small compared with t. This allows for an equilibration between 
a few traps at every stage. From the asymptotic solutions and definition 
we have 

l ~N~r -h- dh~ t U=0 " Jo dqPd{q)q (M) 

ll<...<lr 11. ..l r 

with Pd{q') tending to (cf. eqs.© and flHf )) 

P d (q')=x5(q') + (l-x)5(q'-q). (15) 

Assuming a Boltzmann distribution restricted to the visited region of phase 
space at any stage, and that the long-time traps tend to verify clustering, 
starting from eqs.(|l^-(|T5"l) we can use similar arguments to those leading to 
the interpretation of the static P(q) [Q] in order to understand the dynamical 
Pd{q) (note that in this model Pd{q) 7^ P{q))- Much of the interpretation of 
the static P(q) carries on to the measure associated to the pseudo-equilibria. 
There is however an important difference: the identity of the dominating 
pseudo-states changes with time (otherwise these would be bona fide states 
which we have previously seen they are not). The dynamic phase transition 
takes place when x reaches one in a manner that resembles the static transi- 
tion P; at that point the threshold energy coincides with the paramagnetic 
energy. 

It can come as a surprise that analytic results can be obtained at all in 
such non-equilibrium situations: the underlying reason is the weakness of the 
memory of the system. It would be interesting to understand whether this 
also holds for more realistic systems. 

We wish to acknowledge useful discussions with A. Crisanti, and H. 
Rieger. We are indebted with S. Franz, E. Marinari, G. Parisi and M.A. 
Virasoro for critical reading of the manuscript and suggestions. 
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